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Abstract 

In this paper we show that a universe dominated by two components, namely domain walls and 
dissipative dark matter such that the former slowly decays into the latter may drive power law 
cosmological acceleration consistent with the high redshift supernovae data, while the ratio between 
the energy densities of both components remains fixed at late times thus solving the coincidence 
problem of present acceleration. Likewise, we estimate the aforesaid ratio at early times. 
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I. INTRODUCTION 



In the last years measurements of the luminosity-redshift of type la supernovae have 
provided growing evidence that our Universe has recently entered a phase of accelerated 
expansion 0] . In Einstein gravity accelerated expansion of a spatially flat universe requires 
that the cosmological dynamics be dominated by some exotic component with a negative 
pressure high enough to violate the strong energy condition. 

The position of the first acoustic peak of the CMB as witnessed by the WMAP satellite 

combined with large scale observations of mass distribution suggests that this may be 
the case as the spatial geometry is nearly flat while the mass density in clustered matter is 
roughly half the critical density PI . This is also consistent with studies of the baryon fraction 
in galaxy clusters [5] together with big bang nucleosynthesis calculations 0]. Therefore, some 
additional non-clustered form of energy (dark energy) is needed to attain critical density. 

The precise nature of the dark energy component has become a subject of intense de- 
bate [t^. Dark energy candidates include the cosmological constant (or vacuum energy), 
quintessence (a nearly spatially homogeneous but time-dependent scalar field) 0, tachyon 
field joL Chaplygin gas [10], solid dark energy ^|, a dense network of low-tension domain 
walls 13; among others. 

In introducing dark energy as a novel component most of these candidates face the em- 
barrasing question: "Why are the energy densities of both components of precisely the same 
order today?" This is the so-called "coincidence problem" jl^ . 

The aim is this paper is to study a model of late acceleration dominated by two compo- 
nents, namely a strongly frustrated network of domain walls and dissipative dark matter. 
It is usually assumed that domain walls and dark matter evolve independently. An inter- 
esting alternative is to consider that the evolution of the domain walls is slowed down by 
their interaction with the dark matter fluid ^|. In that model the average wall velocity is 
determined by the ratio of wall and dark matter energy densities. In this respect, it seems 
worthwhile to explore a model featuring a coupling between both components. Here, we 
assume that the two components do not conserve separately but the domain walls slowly 
decay in dark matter. This process could be realized by a broken symmetry that gives do- 
main walls and then the symmetry could be restored and the walls would slowly decay . 
Here, by slowly we mean that pdw/ Pdw < H, where pdw is the energy density associated to 
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the domain walls (subscript dw), H = a/a is the Hubble factor, a{t) is the scale factor of 
the Friedmann-Robert son- Walker metric and the over dots denote derivative with respect 
to cosmic time. In this case the energy in the walls would have to go into dark matter 
particles. As it turns out, this scenario is compatible with a stage of accelerated expansion 
governed by a near power law scale factor consistent with the high redshift supernovae data. 
In addition, the ratio between the energy densities pm/ Pdw stays constant at late times. This 
solves the coincidence problem that afflicts many approaches to late acceleration without 
recourse to fine tuning. 

Domain walls are ubiquitous in field theory models with spontaneously broken discrete 
symmetries ITj] and they appear to be compatible with the observations of large scale 
structure 18[, and the CMB anisotropy spectrum provided that the domain walls network 
is strongly frustrated In such a case its equation of state parameter w = Pdw/ pdw may 
vary with expansion according to 

2 a din rj 

3 3 (iln a' 



where rj the conformal time, dt = a{t)dri and a{< 1) a constant 13|. Further, Conversi et 
al. have recently shown that even if the network is not frustrated, no conflict with the 
CMB anisotropics arises provided the reduced Hubble parameter is less than 0.65 (as some 
experiments seem to indicate [2^) and the matter density parameter lies above 0.35 (as the 
case might be, see Ref. 3])- In section II we explore the possibility that domain walls may 
act as dark energy in late expansion in such a way that the coincidence problem may be 
solved. This analysis extends previous studies of two interacting components -quintessence 
and dark matter-, in whichja^constant stable ratio arises between them at late times due 
to their mutual interaction [22], to a universe dominated by domain walls and dark matter. 
Obviouly, the ratio pm/ Pdm cannot have stayed constant during the whole Universe history, 
otherwise the nucleosynthesis, CMB and cosmic structure scenarios would have been very 
different. This is why section HI focus on the evolution of that ratio and evaluates it at the 
time of domain walls formation. Section V calculates the stateflnder parameters. Finally, 
section V summarizes our flndings. 
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II. DOMAIN WALLS-DARK MATTER INTERACTION 



We consider a spatially flat Friedmann-Lemaitre-Robertson-Walker universe dominated 
by two interacting components, a strongly frustrated network of domain walls and dark 
matter (subscript m). In the following we shall assume that the domain walls behave as a 
fluid as a whole. Consequently, the stress-energy tensor of the cosmic medium takes up the 
hydrodynamic form and the corresponding Friedmann equation is 



K 



= 8nG), 



(1) 



where p = Pm + Pdw is the total energy density. 

Since by assumption both components do not conserve separately a source/loss term 
(say, Q) must enter the energy balances 



Q, 



(2) 



and 



Pdw + ^H{pdw + Pdw) = -Q, 



(3) 



where tt^ stands for the dissipative pressure in the dark matter fluid. Bar ring superfluids, 
the dissipative pressure arises rather naturally in all fluids in Nature I23I . 24 1 whereby we 
have no a priori reason to ignore it (unless we knew that dark matter were in a superfluid 
state, something that does not seem to be the case). Likewise, a comparatively high negative 
pressure may arise within the dark matter fluid by the decay of dark matter particles into 
dark particles j^]. The decay of dark particles in the halo of galaxies has been considered 
yy different authors to explain the absence of steep cusp in the halos of the galaxies -see 
26[ and references therein. Further, the development of density inhomogeneities has been 
pointed out as a source of negative pressure in the cosmic fluid |27|. 
We are interested in studying the dynamics of the ratio Pm/ Pdw 
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By virtue of equations ^ and Q we can write 



d_ 

dt 



Pm 


= 3H^ I 


_Pdw_ 


Pdw \ 
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Pm Pdw "iH p„ 



(5) 



where for simplicity we assumed that the matter fluid behaves as pressureless dust. 
A stationary stage, i.e., d{pm/ Pdw) / dt = is achieved for 



w 



Q = —^H— Pdw 

1 + ro 



Tlr. 



W Pdu 



(6) 



where ro stands for the stationary value of the ratio pm/ pdw at present time. 

To study the stability of this stationary solution against small perturbations we introduce 



the ansatz — = ro + e into Eq.® to get 

Pdw 



e = 3H 



w + 



P Q_ _TWn 
pdw pm 

3H 

Pm 



ro + e 



(7) 



Following [28[ we assume that the interaction and the dissipative pressure are related to the 
total energy density hj Q = Sc^Hp, and tt^ = —b'^p, where and 6^ are positive-definite 
constants. (It should be borne in mind that vr^ is always negative for expanding fluids 
Thus, up to first order in e one obtains 



3H 



c\rl - 1) - 
ro 



(8) 



Therefore the stationary solution will be stable for 



ro < \/ 1 + 



62 



(9) 



which is compatible with the observation that pm/ Pdm < 1 at present. 
From (jH)) it follows that 



ro 



(1 + ro)^ 



w + h 



2 1+^0 
ro 



(10) 
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This sets on the dissipative parameter the upper bound 6^ < — rou;/(l + tq) since as noted 
above must be positive. 

From Eqs.Q, (jSI) and (jH)) along with the assumption that w varies slowly (or that it is 
piecewise constant) the dependence of the energy densities on the scale factor at late times 
is readily found 

p-m OC a"^, pdu, oc a"^, (11) 



where 



/3 



w 



l+ro 



+ 1-6^ 



(12) 



Combining Eq. (fTT|) with Friedmann's equation we get 



a (X t 



2/0 



(13) 



whereby accelerated expansion will occur for (3 < 2. 

When the Universe is dominated by just one component with constant equation of state 
parameter w^:, the scale factor varies as 

^2/3{i+«,.)_ Therefore, in our case, the effective 

parameter is 

l + ro 

For the sake of numerical evaluation let us consider tq = 3/7 -see e.g. Ref. For 
w = —0.5 and b"^ = 0.14, one follows We// = —0.49, whence /3 = 1.53; while for w = 
—0.66 and 6^ = 0.19, one obtains Wejf = —0.66, whence /3 = 1.03. These Weff figures are 
only marginally compatible with analyses that tend to push We// closer and closer to the 
cosmological constant value and even beyond it 0,13 • However, it should be noted that, 
as observed by Maor et al. [30], these analysis are based on combining luminosity-redshift 
measurements with models that consider the equation of state parameter of the dark energy 
Wx constant which automatically introduces a bias that favors larger negatives values of 



26 




0.0 0.2 0.4 0.6 0.8 1.0 

Redshift 



FIG. 1: Effective magnitude vs redshift for the model discussed in the paper. To plot the graphs 
the expression m'^J^ = M^f + 5log{HQdL), with Msf = —3.9 and d/, = (1 + z) H^^{z')dz' , was 
used. For comparison, we include the prediction of the ACDM model with pm = 0.3. The data 
were borrowed from figure 6 of Ref. j^. 

than if it varied with expansion. On the other hand, as Fig. shows the agreement of our 
model with the supernovae data of Knop et al. is excellent. 

At this stage, one may wonder why at all one needs a coupling between both components 
and a dissipative bulk pressure in the dark matter if the the most negative value for We/f 
essentially coincides with the one provided by domain walls (—2/3). The answer is that 
these ingredients are key to guarantee that the ratio pm/ Pdw will remain constant at late 
times. In our view, models of late acceleration that fail to address the coincidence problem 
are not very useful. 
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III. THE RATIO prn/pdw AT EARLY TIMES 



Obviously, at the time of domain walls formation (about 10~^^ seconds after the big 
bang j^l) the ratio between the energy densities must have been very different from its 
present value, tq. In this connection it is mandatory to study the evolution of r = Pm/ Pdm 
with the Universe expansion if one wishes to evaluate r at that epoch. Here we shall follow 



a procedure parallel to that of Ref. 



The evolution equation for r is given by Eq.(|H)), 



with Q = 3c^Hp and TCr, 



-b'^p. Its stationary points. 



b"^ + w 
2c2 



(14) 



where 



^ w"^ w f b"^ \ b'^ 



follow from setting the right hand side of Eq.© to zero with w = constant. Notice that 11 
determines the difference between the stationary solutions r"*" — r~ = 2-\/TT, and that their 
sum is r"^ + = —2 — {b^ + w)c~'^ . To these relations the constraint 

r+r' = 1 + ^ (16) 



is to be added. The stability of the stationary points is determined by the sign of [dr /dr\r±. 
It can be seen that that is unstable and r~ is asymptotically stable. This tells us that the 
solution of Eq.(jS)) evolves from to r~ with > r~ -see Fig.Q In order to describe the 
transition from the decelerated phase to the accelerated one we resort to the deceleration 
parameter q = —aa/a?. Before domain walls dominance the Universe is decelerating, g > 0, 
and later changes to g < 0. In virtue of the Friedmann equation and the evolution 
equation for the Hubble factor 

H = --[pm+ {.^ + W)pdu, + T^m], (17) 
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FIG. 2: The ratio between energy densities Pm/pdm decreases with expansion up to an asymptotic 
value at the present epoch. In drawing both curves we assumed = 0.14 and = 0.04. 



we can write 



^=2 



1 + 



w 



1 + r 



(18) 



It is worthy of note that the condition g+ = > > q{r ) — q , imphes the constraint 



w 1 y2 w 1 



(19) 



Likewise, it is convenient, for later purpose, to estimate (? in terms of and r . From 
Eqs.ffT^. ([T^ and the constraint (fT?)|) . we readily obtain 



1 - Iq- 



3r+(l + r 



(20) 



in agreement with the previous statement that (? ought to be a positive-definite quantity. 



The condition q{rac) = yields the density ratio at the beginning of the acceleration 



-3w + 362-1 



(21) 



Using the restriction 36^ < 1 we get > — (1 + Sw). Obviously r+ > r^c > r~, as can be 
checked. 

The observational constraints of type la (SNela) supernovae may be used to determine 



le corrected 
l|. We may 



the values of and b'^. For each supernova, its redshift Zi, dispersion o"j and t 
magnitude rrii, were computed, resulting a ~ as the most likely value of a 
tentatively assume that for cr = the Universe is settled in an asymptotic energy density 
ratio r~ ~ tq, with constant deceleration parameter . On the order hand, for power-law 
expansion given by Eq. (fT3|) . we find rf^ oc (1 + z)[{l + zY — 1], with 5 = —q~. Thus, we get 
< (5 < 1, its most likely value being 6 ~ 0.395. 

Let us now study the implications of the above results for and b"^. Using Eq.(|20p and the 
constraint (fT^ we obtain 



c' + b' 



1 + r- 



2q- 



1 + ro 



1 + 2(5 



0.179. 



(22) 



To assign a value to r+ (which may be identified with the value of r at the end of domain 
walls formation), we assume it should not be far from the r value at the epoch of primordial 
nucleosynthesis, r^- Note that the nucleosynthesis epoch is the earliest one on which we 
have reliable quantitative information and that the time elapsed between domain walls 
formation and nucleosynthesis is larger by many magnitude orders than the time elapsed 
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between the latter and the present epoch. As can be checked (see the Appendix of Ref . |28l| ) . 



r+ ^ > ^ 6.143, (23) 

' '-^ max 



where the maximum variation in the effective number of neutrino species satisfies 
A(7*(T/v) < 7ANmax/4:, with ANmax = 1 and g*{Tj^) = 10.75 for the standard model 
with three massless neutrinos 0|. Assuming r+ ^ r^v and ~ tq = 3/7 in Eq. ()20j) . it 
follows that < 0.068. Accordingly, 6^ must lie in the range 0.111< < 0.179, whereby 
< < 0.068. From the constraint (|16|) we obtain 

^~ r^ro-1 > 1.633, (24) 

suggesting that the viscous effects are rather important. Note that because the time of 
domain walls formation precedes primordial nucleosynthesis and r decreases with expansion, 
r^v just represents a lower bound on the value of r at that time. 

IV. STATEFINDER PARAMETERS 

Due to the ample choice of dark energy models the deceleration parameter is unable 
by itself to discriminate between them nor even from the ACDM model. This is why 
Sahni et al. introduced two dimensionless, geometrical parameters, the so-called 

statefinder pair (r, s), to more fully characterize the different models and in particular to 
help discriminate between the ACDM model and those based on a varying equation of 
state. These parameters are defined as 



a _ r — 1 

o^' ^ " 3(q- 



1. ■ (25) 

2) 



In our case, we find that at late times 
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being related by r = 1 + — 1). 

In more general situations such that w may not be reasonably well approximated by a 
constant, the deceleration parameter and the statefinder pair can be calculated from the 
Friedmann equation and the evolution equation for the Hubble factor given by Eq. lfTTj) . 
The result is now 

g = i [1 + ^{wn,^ - 62)] , (27) 
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dw 



(1 _ ^2) [^2(^ + 1) + 1] _ (1 + + ^) _ 



(28) 



+ 



rw 
r + 1 



w 



w 
3H 



and 



dw 



(62 -wildu,) 



I (1 - b^)[b\r + 1) + 1] - (1 + ^)[(1 + w)- b^] (29) 



b'^w rw"^ w 

+ + 7+^+7777 

r r + 1 SH 

where Qdw = Pdw/ {"^H"^) is the density parameter of the domain walls component. 



V. DISCUSSION 



We have presented a cosmological model of late acceleration driven by domain walls 
coupled with dissipative dark matter fully consistent with the high redshift supernovae data. 
At late times the coupling leads to a stable constant ratio between the energy densities 
of both components compatible with accelerated expansion, thus solving the coincidence 
problem without recourse to fine tuning. This ratio can be thought as constant at late times 
only, thereby we have studied its evolution with the Universe expansion and estimated 
its value, r^v, at the nucleosynthesis epoch. This sets a lower bound on the value of r 
at the time of domain walls formation. In so doing we have implicitly assumed that the 
interaction domain walls-dark matter starts right after the former component comes into 
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existence. In a more detailed analysis the radiation component should enter the picture 
as well as it dominates the expansion from the end of the inflationary era till matter- 
radiation equality. However, this phase is extremely short compared with the matter plus 
dark energy dominance period. We believe its influence should not substantially alter the 
main conclusions of the paper. 

We have, moreover, calculated the statefinder parameters of this model. It is to be hoped 
that in some not distant future observational techniques able to measure the (r, s) pair will 
be in place so that these two parameters will become useful tools to ascertain the nature of 
dark energy. 
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